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ABSTRACT 

I present a catalogue of positions and correlated flux densities of 109 compact 
extragalactic radio sources in the Galactic plane determined from analysis of a 48 
hour VLBI experiment at 22 GHz with the European VLBI Network. The median 
position uncertainty is 9 mas. The correlated flux densities of detected sources are in 
the range of 20 to 300 mJy. In addition to target sources, nine water masers have been 
detected, two of them new. I derived position of masers with accuracies 30 to 200 mas 
and determined velocities of maser components and their correlated flux densities. The 
catalogue and supporting material is available at http://astrogeo.org/egaps. 
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1 INTRODUCTION 

The differential very long baseline (VLBI) astrometry be- 
came a powerful instrument for a study of the 3D structure 
and dynamics of our Galaxy. Analysis of differential obser- 
vations allows us to determine parallaxes of compact radio 
sources with accuracies down to several tens of microarc- 
seconds, positions with accuracies 0.05 mas and proper mo- 
tions. The differential astrometry observations are made in 
the so-called phase referencing mode when radio telescopes 
of a VLBI network rapidly switch from a target source to 
a near-by calibrator. Analysis of differential fringe phases 
allows to dilute the contribution of mismodeled path delays 
in the atmosphere by the factor of target to calibrator sep- 
aration in radians, i.e. 20-50 times. Nowadays, observations 
in the phase-referencing mode are widely used for imaging 
faint sources. The atmosphere path delay fluctuations limit 
the integration time, depending on frequency, to 0.2-10 min- 
utes. Phase referencing observations allow to overcome this 
limit and to integrate the signal for hours. According to 
IWrobell (|2009l ), 63% of Very Long Baseline Array (VLB A) 
observations in 2003-2008 were made in the phase referenc- 
ing mode. 

However, the feasibility of phase-referencing observa- 
tions is determined by the availability of a pool of cal- 
ibrators with precisely known positions. The first cata- 
logue of source coordinates determined with VLBI con- 
tained 35 objects l|Cohen fc Shaffer~lll97ll ). Since then, hun- 
dreds of sources have been observed under geodesy and 
astrometry VLBI observing programs at 8.6 and 2.3 GHz 
(X and S bands) using the Mark3 recording system at 
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the International VLBI Service for Geodesy and Astrom- 
etry (IVS) network. Analysis of observations in 1980s 
and 1990s result ed in the ICRF catalogue of 608 sources 
(|Ma et al.l Il998l ). Later, using the VLBA, positions of 
over 5000 other compact radio sources were determined 
in the framew o rk of t he VLBA Calibrator S urvey (VCS) 
ilBeaslev et alj 120021: iFomalont et all 120031: JPetrov et all 
l2005l . 120061 ; iKovalev et al.1 120071; IPetrov et alj|2007al). the 
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VIPS program (|Helmboldt et al.l _ 

l2011dl ). the BeSSel Calibrator Survey (|lmmer et al.ll2011n T 
and in the number o f on-going program s: the geodetic 
VLBA program RDV l|Petrov et al.ll2009l ), the program of 
a study of active galaxy nuclea (AGNs) at parsec scales de- 
tected with Fermi (Kovalev & Petrov, paper in prepara- 
tion), the program of observing radio- loud 2MASS galax- 
ies (Petrov, paper in preparation ), and the program of ob- 
serving optic ally bright quasars l|Bourda et al.l 120081 . 1201 ll : 
|Petrovll201ld ). The use of the Australian Long Baseline Ar- 
ray (LBA) ex tended the catalogue of calibrators to the zone 
[-90°, -40°] (IPetrov et al.ll2011bl ). 

By April 2011, the total number of radio sources with 
positions determined with the VLBI in the absolute astrom- 
etry mode reached 6123 and it continues to grow. On aver- 
age, the probability of finding a calibrator within a 3° radius 
of a given position is 97%. However, the density of calibra- 
tors in the Galactic plane is still low and an increase of the 
number of calibrators is badly needed for many VLBI Galac- 
tic astronomy projects. Finding extragalactic sources visible 
through the Galactic plane region is more difficult for sev- 
eral reasons. First, the region is crowded with many galactic 
objects. Combining independent low-resolution observations 
at different frequencies makes determination of their spec- 
tra problematic due to a risk of source mis-identification. 
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Second, many potential candidates with flat spectra are ex- 
tended galactic objects, such as planetary nebulae or com- 
pact HII regions, that cannot be detected with VLBI. Third, 
the apparent angular size of extragalactic objects observed 
through high plasma density near the Galactic plane is 
broadened by Galactic scattering and cannot be detected 
at low frequencies on baselines longer than several thousand 
kilometers. 

To address the problem of increasing the density of cal- 
ibrators within the Galactic plane, we made in 2005-2006 
a dedic ated blind fringe su rvey of 2496 objects with the 
VERA (|Petrov et alj|2007bh radio interferometer at the K- 
band (22 GHz). Those sources which were detected with the 
VERA were re-observed with the VLBA in 20 06 at K-band 
in the framework of the VGaPS campaign (|Petrov et al.l 
l2011al ). These observations allowed to determine positions 
of 176 new sources within 6°of the Galactic plane with the 
median accuracy of 0.9 mas. The detection limit of VERA 
was 200-300 mjy, the detection limit of the follow-up VLBA 
observations was 70-90 mjy depending on a baselin^E In- 
dependently, a team led by Mark Reid searched for Galactic 
plane calibrators in 2010 using the VLBA in the framework 
of BeSSeL project and reported detection of 198 sources 
llmmer et al.l (|20 1 lT ) . 82 of them are new. 

In this paper I present results from a 48 hour experiment 
at the European VLBI Network (EVN) observed in October 
2009 called the EVN Galactic Plane Survey (EGaPS). The 
goal of this experiment is to further increase the calibrator 
source density in the zone within 6° of the Galactic plane and 
with declinations > —20°. Using highly sensitive antennas, 
I was able to detect sources as weak as 20 mjy, a factor of 4 
weaker than in the prior VGaPS campaign. These calibra- 
tors are aimed to be used for Galactic astronomy projects. 
The selection of candidate sources and the scheduling strat- 
egy is discussed in section The station setups during the 
observing sessions is described in section [3] The correlation 
and post-correlation analysis is discussed in section [4] The 
catalogue of source positions and correlated flux densities 
is presented in section [6] An error analysis of observations, 
including evaluation of systematic errors, is given in subsec- 
tion 231 an d the results are summarized in section [7] 



2 CANDIDATE SOURCE SELECTION 

In the past, the list of bright compact sources detected at the 
Very Large Array (VLA) and/or MERLIN served as a pool 
for candidates for VLBI calibrator surveys. Almost all these 
objects have been already surveyed with the VLBA. Since 
the combined catalogue of sources observed in the absolute 
astrometry mode is complete at the 150-200 mjy level, new 
calibrators should be relatively faint objects. Most all-sky 
surveys at frequencies higher than 2 GHz are either incom- 
plete and/or not deep enough. As a result, information about 
spectral indices is sparse and is often unreliable. For this rea- 
son, it is more difficult to find remaining good flat spectrum 
candidate sources. 



1 The K-band receivers at the VLBA stations were upgraded in 
2007, after these observations, which improved the sensitivity of 
the array at 24 GHz at the sam e recording rate by a factor of 2 
according to lWalker et alj j200ST) . 



Figure 1. The EVN stations that participated in the EGaPS 
experiment. The longest meridional baselines NOTO/ONSALA60 is 
2280 km long. The longest longitudinal baseline noto/yebes40m 
is 1616 km long. 
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New calibrator candidates were found by analyz- 
ing 395 radio astronomical catal ogues using the CATS 
databa se dVerkhodanov et al.|[l997h . T he database include s 
NVSS dCondon et al. Ill998h. CLASS (|Mvers et al. 1 120031 ') . 
CRATES (iHealev et al l l2007t ), ATCA survey at 20 GH: 



jMurphv et alj|2010h . and many other catalogues. Initially, 
I selected 1100 sources which satisfy the following criteria: 

• were not observed previously with VLBI; 

• have galactic latitude |6| < 6°; 

• have declination > —20°; 

• have at least two measurements of their flux density 
that allow estimation of their spectral indices at frequencies 
higher than 2 GHz; 

• have single-dish flux densities extrapolated to 22 GHz 
> 80mJy; 

• have spectral indices a flatter than —0.5 (S ~ v a ). 

The majority of these sources w ere selected by cross - 
identification of the NVSS and GB6 l|Gregorv et al. ||l996h 
catalogues. Then I scrutinized the source spectra and re- 
jected approximately 50% of the sources with unreliable 
spectra or the objects known as galactic sources. The re- 
maining list contains 559 candidates. 

2.1 Observation scheduling 

The following stations were scheduled for 48 hour obser- 
vations: EFLSBERG, JODRELL2, MEDICINA, NOTO, ONSALA60, 
and yebes40m (see Figure [T] The observation schedule was 
prepared with the software program sur_sked. The schedul- 
ing goal was to observe each target source at all antennas of 
the array in 2 scans of 120 seconds each. 

The scheduling algorithm found the sequence of sources 
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that minimized slewing time. The minimum time between 
consecutive observations of the same source was set to 
4 hours. For deciding which source to put into the sched- 
ule, the algorithm calculated the elevation of each candidate 
object, the slewing time, the time interval between previous 
observation of the same source, and the ratio of the remain- 
ing visibility time to the total visibility time. For each source 
with the elevation exceeding 15°above the horizon, a score 
was computed. The score is the function of slewing time, the 
ratio of the remaining visibility time to the total visibility 
time, and the time interval from the previous observation. 
A source with the highest score was put into the schedule 
and then the procedure was repeated. 

Every 1.5 hr a burst of 4 stron g compact sources with 
known maps from the K/Q survey riLanvi et al.ll2O10D was 
observed: two strong objects at elevation angles 10-30° and 
two strong objects at elevations 50-90°. The purpose of 
scheduling these calibrators was a) to allow estimation of 
the troposphere path delay in zenith direction; b) to eval- 
uate the atmosphere opacity; c) to use them for complex 
bandpass calibration; d) to provide a strong connection be- 
tween the new catalogue to the old catalogue of compact 
sources. 

The schedule had 369 target objects and 69 calibrators. 
Among target sources, 344 objects were observed in 2 scans 
and 25 sources were observed in 1 scan. 



3 OBSERVATIONS 

Observations took place on 2009 October 27-29. Eight in- 
termediate frequencies (IF), both upper and lower sub-band 
in the range of [22.09999, 22.35599] GHz were observed in 
a single left circular polarization with two bits per sample. 
The data were recorded with the Mark-5 system with the ag- 
gregate rate 1024 Mbit/sec. The contiguous frequency setup 
is far from optimal for an absolute astrometry experiment 
that utilizes group delays, since the uncertainties of a group 
delay at a given signal-to-noise ratio (SNR) is reciprocal to 
the root mean squares (rms) of IFs. Two stations, medicina 
and ONSALA60 could record within 720 MHz, three stations 
equipped with the VLBA data acquisition system, eflsberg 
and NOTO yebes40m, could record within 500 MHz and JO- 
drell2 could record within 160 MHz. Spreading IFs over 
500 MHz would reduce group delay uncertainty by a factor 
of 2.4 at all baselines, except jodrell2 where it would in- 
crease the uncertainty by a factor of 2 since some channels 
will be missing. However, in order to record a band wider 
than 256 MHz, the local oscillator (LO) frequency should be 
changed with respect to those that are used in other EVN 
experiments. This would require manual intervention. When 
I realized that for logistical reasons there is a high risk that 
this change would not be made at some stations, which could 
ruin the experiment, I fell back to the contiguous frequency 
setup. 

Station yebes40m lost first 7 hours of the experiment 
because of the Mark-5B failure. The system temperature 
of YEBES40M was around 210 K during next three hours 
because the vertex was closed. It dropped to the normal 
range of 60-100 K when the operator noticed it and opened 
the vertex. Station JODRELL2 did not show fringes in IF 11- 
IF16. 



4 DATA ANALYSIS 

Analysis of the EGaPS data is similar to that performed 
for the VLBA Galactic plane calibrator survey. Detailed 
description o f the analysis technique can be found in 
iPetrov et ail (|2011aT ). Here the analysis procedure is briefly 
outlined. 



4.1 Data correlation and post-correlation analysis 

The data were correlated at the Bonn DiFX software corre- 
lator iDeller et alj l|2007h with the accumulation periods of 
0.12 s and with the spectral resolution of 125 KHz. These 
correlation parameters provided the field of view over 1', 
comparable with the width of the beam of eflsberg radio 
telescope. I anticipated the a priori source coordinates of 
some sources may be wrong at the arc-minute level due to 
source mis-identification. 

The correlator generates the spectrum of the cross- 
correlation function that was written in the format c ompli- 
ant with the FITS-IDI specifications (|GreiserJ|2009T ). Data 
analysis was performed with software VTMJvQ At the first 
step, the fringe amplitudes were corrected for the signal dis- 
tortion in the sampler. Observations of strong sources were 
used for deriving complex bandpasses. Since in 2010 the 
DiFX correlator did not have an ability to extract phase 
calibration signal from the data, no phase calibration was 
applied. After that, the group delay, phase delay rate, group 
delay rate, and fringe phase were determined for each scan 
at each baseline using the fringe fitting procedure. These es- 
timates maximize the sum of the cross-correlation spectrum 
coherently averaged over all accumulation periods of a scan 
and over all frequency channels in all IFs. After the first run 
of fringe fitting, 8 observations at each baseline with the 
strongest SNR were used to adjust the station-based com- 
plex bandpass corrections, and the procedure of computing 
group delays was repeated with the bandpass applied. 

Then the results of fringe fitting were exported to 
the VTD/post-Solve VLBI analysis softwar^fjfor interactive 
processing group delays with the SNR high enough to ensure 
the probability of false detection less than 0.001. Analysis 
of the probability distribution function of achieved SNR al- 
lowed me to find that the probability of a false detection in 
our experiment is less than 0.001 when SNR> 6.14. Theoret- 
ical path delays were computed according to the state-of-the 
art parametric model as well as their partial derivatives. The 
model included the contribution of the ionosphere to group 
delay and phase delay rate computed using the total electron 
contents maps derived from analysis of Global Positioning 
System (GPS ) observations by the analysis center CODE 
l|Schaerlll998l l. Small differences between group delays and 
theoretical path delay as well as between the measured and 
theoretical delay rate were used for an interactive estima- 
tion of corrections to a parametric model that describes the 
observations with least squares (LSQ). Coordinates of tar- 
get sources, positions of all stations, except the reference 
one, parameters of the splines that describe corrections to 
the a priori path delay in the neutral atmosphere in the 
zenith direction for all stations, and parameters of another 
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Figure 2. The distribution of calibrator sources within 10° of the Galactic plane. The filled disks denote sources detected in the EGaPS 
campaign. The hollow circles denote calibrator sources known from previous observations. New sources from the BeSSeL VLBI survey 
are not shown. 



=+60° 6=+40° 6=+20° 6=0° 




180 210 240" 
Galactic longitude (deg) 



splines that describe the clock function were solved for in a 
single LSQ solution using group delays. Outliers were iden- 
tified and temporarily suppressed, and additive corrections 
to weights of observables were determined. Then the fringe 
fitting procedure was repeated for outliers with phase de- 
lay rates and group delays evaluated in a narrow window 
around the expected value computed on the basis of results 
of the previous LSQ solution. New estimates of group delays 
for points with probabilities of false detection less than 0.1, 
which corresponds to the SNR> 4.9 for the narrow fringe 
search window, were used in the next step of the interac- 
tive analysis procedure. For observations detected with the 
narrow window search the status outlier was lifted and they 
were selected back for using in further analysis. Parame- 
ter estimation, elimination of remaining outliers and adjust- 
ments of additive weight corrections was then repeated. In 
total, group delays of 2471 observations out of 11675 sched- 
uled were used in the solution. There were 178 sources de- 
tected in 3 to 28 observations, including 109 target sources 
and 69 calibrators. 

4.2 Source position determination 

Results of the interactive solution provided a clean dataset 
of group delays with updated weights from the EGaPS ex- 
periment. The dataset that was used for the final parameter 
estimation utilized all dual-band S/X data acquired under 
absolute astrometry and space geodesy programs from April 
1980 through February 2011 and included the K-band data 
from the EGaPS experiment, in total 8 million observations. 
The estimated parameters are right ascensions and declina- 
tion of all sources, coordinates and velocities of all stations, 
coefficients of B-spline expansion of non-linear motion for 
16 stations, coefficients of harmonic site position variations 
of 48 stations at four frequencies: annual, semi-annual, diur- 
nal, semi-diurnal, and axis offsets for 67 stations. They were 
adjusted using all the data. Estimated parameters also in- 
cluded Earth orientation parameters for each observing ses- 
sion, parameters of clock function and residual atmosphere 
path delays in the zenith direction modeled with the linear 
B-spline with interval 60 and 20 minutes respectively. All 
parameters were estimated in a single LSQ run. 

The system of LSQ equations has an incomplete rank 
and defines a family of solutions. In order to pick a spe- 
cific element from this family of solutions, I applied the 
no-net rotation constraints on the position s of 212 sources 
marked as "defining" in the ICRF catalogue (|Ma et al.lll99&t ) 
that required the positions of these source in the new cat- 
alogue to have no rotation with respect to the position in 
the ICRF catalogue. No-net rotation and no-net-translation 
constraints on site positions and linear velocities were also 
applied. The specific choice of identifying constraints was 



Table 1. Estimates of JODRELL2 positions in meters on epoch 
2009.10.27 

X 3822846.633 ± 0.026 
Y -153802.071 ± 0.009 
Z 5086286.064 ± 0.036 

made to preserve the continuity of the new catalogue with 
other VLBI solutions made during last 15 years. 

The global solution sets the orientation of the array with 
respect to the ensemble of > 5000 extragalactic remote radio 
sources. The orientation is defined by the continuous series 
of Earth orientation parameters and parameters of the em- 
pirical model of site position variations over 30 years eval- 
uated together with source coordinates. Common sources 
observed in the EGaPS experiment as amplitude calibrators 
provided a connection between the new catalogue and the 
old catalogue of compact sources. 

A valuable by-product of analysis of EGaPS observa- 
tions is estimates of jodrell2 positions (see Table [TJ. This 
is the second experiment with participation of this station 
under astrometry or geodesy programs. The previous EVN 
experiment TP001 carried out on 2000 November 23 was 
made at 5 GHz with the bandwidth spread over 108 MHz 
l|Charlot et alj [2002). Comparison of results of reanalysis 
of TP001 and EGaPS made using exactly the same a pri- 
ori model and parameter estimation technique showed that 
the differences between jodrell2 position reduced to the 
J2000.0 epoch using the velocities predicted by the N NR- 
NUVEL1A plate tectonic model (|DeMets et all I1994T I are 
38, 20, and 9 mm at X, Y and Z coordinate respectively. Al- 
though position of jodrell2 from the 5 GHz observations 
with formal uncertainties 10-20 mm suffered from residual 
errors in modeling path delay through the ionosphere which 
are 20 times greater than at 22 GHz, the differences in po- 
sitions are within 1-2 reported uncertainties. 

4.3 Position error analysis 

The formal uncertainties of semi-major axis of the error el- 
lipses of position estimates of 109 target sources range from 
2 to 60 mas with the median value of 8 mas. They are based 
on the error propagation law of group delay uncertainties 
derived by the fringe fitting procedure. 

Including in observing schedule a considerable number 
of calibrator sources with positions known at the sub-mas 
level facilitated the error analysis. In order to evaluate the 
level of systematic errors, I used a similar technique that 
was developed for processing VGaPS observations. The list 
69 calibrators was ordered in increasing their right ascen- 
sions and split it into two subsets with even and odd in- 
dices. I made two special global solutions using all the data 
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from 1980 through 2011. I excluded in solution A calibra- 
tors with even indices from all experiments, but the EGaPS. 
I excluded in solution B calibrators with odd indices from all 
experiments, but the EGaPS. I compared positions of cal- 
ibrator sources from solutions A and B with their position 
from the dual-band global solution C that used all sources 
in all experiments, except the EGaPS. Positions of sources 
from solution C is determined with accuracies 0.1-0.3 mas 
and can be considered as true for the purpose of this com- 
parison. 

The wrms of the differences of 69 calibrator sources 
from special solutions A and B and the reference solution 
C are 1.8 mas in declinations and 3.7 mas in right ascen- 
sions scaled by cos 5 with \ 2 P er degree of freedom 1.3 and 
3.3 respectively. If to add in quadrature to uncertainties in 
right ascensions and declinations the noise with the standard 
deviations 4.0 m&s/cosS and 1.0 mas respectively, the x 2 per 
degree becomes close to 1. Therefore, I inflated the reported 
uncertainties in right ascensions and declinations of target 
sources by 4.0 m&s/cosS and 1.0 mas in quadrature. 



4.4 Data analysis: correlated flux density 
determination 

Each detected source has from 3 to 30 observations, with 
the median number of 8. The dataset is too sparse to pro- 
duce meaningful images. In this study I limited my analysis 
with mean correlated flux density estimates in two ranges of 
lengths of the baseline projections onto the plane tangential 
to the source without inversion of calibrated visibility data. 
Information about the correlated flux density is needed for 
evaluation of the required integration time when an object 
is used as a phase calibrator. 

Analysis of system temperature measurements revealed 
variations with time and with elevation angle. The measured 
system temperature is considered as a sum of two terms: the 
receiver temperature T r and the contribution of the atmo- 
sphere: 

-P m(e)-| 



Tsys — T r 



Tatm [1 — 1 



(1) 



where T at , n is the average temperature of the atmosphere, 
P is the atmosphere opacity, and m(e) is the wet map- 
ping function: the ratio of the neutral atmosphere non- 
hydrostatic path delay at the elevation e to the atmosphere 
non-hydrostatic path delay in the zenith direction. I omit- 
ted in expression [T] the ground spillover term that was not 
determined for these antennas. Both, the receiver tempera- 
ture and the atmosphere opacity, depend of time. Since in 
astrometric analysis I estimated for each station the non- 
hydrostatic component of the atmosphere path delay in the 
zenith direction which is closely related to the integrated 
column of the water vapor, I can use these estimates for 
modeling time dependence of the opacity. I present the sys- 
tem temperature as 



T sys = T„Bl{t) + T atm [l - e~ {a 



+ i>1"atm(t)) ™(e)l 



(2) 



where a and b are empirical coefficients that relate the opac- 
ity and the estimates of the atmosphere path delay, T atm (t) 
is the estimate of the non-hydrostatic atmosphere path de- 
lay, Bj(t) is the B-spline of the first degree with the piv- 
otal element i, and T r i are the coefficients of the expansion 



the receiver temperature into B-spline basis. I set T atm to 
280K, and evaluated the coefficients of the receiver temper- 
ature and regression parameters a and b by fitting them 
into measurements of system temperatures with the use of 
the non-linear LSQ. Parameter a describes the possible bias 
between estimates of the integrated water vapor contents 
and the non-hydrostatic path delay in the zenith direction. 
The system temperature divided by e"'" + ' ,r »' m ('" m ( e > is free 
from absorption in the atmosphere and equivalent to that 
at the top of the atmosphere. 

The root mean squares (rms) of residuals of the em- 
pirical model of the system temperature are in the range 
of 2-8 K. The time variations of T r were within 10 K at 

EFLSBEREG, JODRELLS2, MEDICINA, and NOTO. The T r was 

around 180 K during first 3 hours at yebes40m and then 
suddenly dropped to 30 K and after that stayed stable. The 
T r at ONSALA60 was unstable during the entire experiment 
and varied in the range 40-250 K. 

The fringe amplitudes were calibrated by multiplying 
them by the system temperature reduced to the top of the 
atmosphere and dividing by the elevation dependent a priori 
gain. 

The a priori antenna gain and/or the term 
r aya /e ( - a+i,Ti ' t ™ (t))m(e) may have a multiplicative er- 
ror. The corrections to antenna gain were evaluated 
by fitting the correlated amplitude to the flux density of 
sources with known brightness distributions. Among sources 
used as amplitude calibrators, brightness distributions are 
publicly availabl^fl for 45 objects from the K/Q survey and 
the VGaPS observing campaigns. 

For each used amplitude calibrator observation with 
known brightness distribution in the form of CLEAN com- 
ponents I predicted the correlated flux density F as 



corr — 



■ (ia x-\-v y) 



(3) 



where c; is the correlated flux density of the ith CLEAN 
component with coordinates x and y with respect to the 
center of the image, and u and v are the projections of the 
baseline vectors on the tangential plane of the source. 

Then I built a system of least square equations for all 
observations of calibrators with known brightness distribu- 
tions used in the astrometric solution: 



Fcorr — yj Qi Qj A C or 



(4) 



and after taking logarithms from left and right hand sides 
solved for corrections to gains g for all stations using LSQ. 
Finally, I applied corrections to gain for observations of all 
other sources. 

The gain correction was in the range 0.7-1.3 for all an- 
tennas, except JODRELL2. The correction to JODRELL2 gain 
was 6.9, and its system equivalent flux density (SEFD) at 
elevations higher than 50°was around 7000 K. 

The detection limit varied significantly between the 
baselines. Sources as weak as 13-16 mjy had the SNR > 6.14 
and therefore, were detected at the most sensitive base- 
line eflsberg/yebes40m. The detection limit at baselines 

4 The database of brightness distributions, correlated flux densi- 
ties, and images of compact radio sources produced with VLBI is 
accessible from http://astrogeo.org/vlbi_images 
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Figure 3. The normalized fringe amplitude of the water maser 
0221+618 at baseline MEDICINA/ONSALA60. The horizontal axis 
shows the frequency offset with respect to the nominal H2O maser 
frequency (6i 6 -5 2 3) 22.235044 GHz. 
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Figure 4. The fringe amplitude of the new water maser 1923+151 
at baseline eflsberg/medioina. The horizontal axis shows the 
frequency offset with respect to the nominal H2 O maser frequency 
22.235044 GHz. 



eflsberg/medicina and eflsberg/noto was in the range 
of 25-30 mjy, at the baseline eflsberg/onsala60 was 
30-40 mjy, and at the baselines eflsberg/jodrell2 and 
medicina/noto was 60-70 mjy. 



5 DATA ANALYSIS: H 2 MASER SOURCES 

During scrutinizing the cross-correlation spectra, I found 
several objects with very strong peaks in the IF9, for exam- 
ple object 0221+618 at baseline medicina/onsala60 (see 
Figure [3]). 

This feature of the fringe spectrum is interpreted as an 
emission from a water maser. The peak frequency of the 
cross-spectrum is shifted within several megahertz due to a 
relative motion of the maser with respect to the geocenter. 
Although this project was not targeted on maser investi- 
gation, nevertheless, I decided to systematically search for 
masers in the data, determine their coordinates and param- 
eters of their spectra. 

Processing the data from narrow-band sources dif- 
fers significantly from processing the data from continuum 
sources. I searched for masers by running 31 trial fringe 
fitting with spectral window of 1 MHz within the IF9 in 
the range of 22.22799-22.24399 GHz. Other spectral con- 
stituents were masked out. The spectral window was shifted 
at 0.5 MHz after each run. Since the maser has emission only 
within a narrow band, suppressing the spectrum beyond the 
search window improved the SNR considerably, because the 
amount of noise is reduced, but the power of the signal re- 
mained the same. I searched for objects with the peak fringe 
amplitude higher than 4 times the spectrum rms. This ap- 
proach helped to discover 9 maser objects. After the sig- 
nal from masers was detected at at least one baseline, the 
spectrum filter was narrowed down even further to reduce 
the contribution of noise away from the spectral lines. This 
allowed me to increase the number of detections at other 
baselines. 

The frequency resolution of the cross-spectrum, 
125 KHz, that correspond to the Doppler velocity resolu- 
tion 1.7 km s _1 , is not sufficient to resolve the water maser 
spectral line. The scans with 7 out of 9 detected masers 




Figure 5. The fringe amplitude of the new water maser 2011+360 
at baseline EFLSBERG/JODRELL2. The horizontal axis shows the 
frequency offset with respect to the nominal H 2 O maser frequency 
22.235044 GHz. 



were re-correlated with 16 times higher spectral resolution of 
7.8 KHz, which corresponds to the Doppler velocity resolu- 
tion 0.1 km s" 1 . Two other masers, 2130+556 and 2247+596 
were discovered after the recorded data were recycled, and 
it was too late to re-correlate them. 

I smoothed the cross-spectrum and fitted the prelim- 
inary model of the amplitude spectrum of maser emission 
A = Ai exp — {(// — fo) 2 /Si } using an iterative non- 
linear LSQ procedure. Here Ai is the peak amplitude of the 
component, fi its frequency, Si is the parameter of the line 
width, and fo is the rest H 2 C> frequency. After determining 
a satisfactory preliminary model of the smoothed spectrum, 
I fitted the final model using the raw cross-spectrum. Since 
masers are not point-like objects, but usually a conglomerate 
of tens or even hundreds spots, the correlated flux density 
at different parts of the uw-plane differs considerably due to 
beatings from different components. Using only 10-20 points 
at the uv-plane is not sufficient for an image reconstruction. 
Therefore, I restricted the amplitude analysis to determina- 
tion of the minimal and maximum correlated flux densities 
of detected components at different baselines. 

The fringe fitting process cannot provide a meaningful 
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estimate of a group delay for a narrow-band source. There- 
fore, I resorted to using phase delay rates for determining the 
source positions. A phase delay rate over 2 minutes of inte- 
gration time at 1024 Gb/s recording rate has the uncertainty 
due to the thermal noise of order 2-10~ 13 /SNR. However, 
unmodeled short-term atmosphere path delay fluctuations 
limit the accuracy of modeling the phase delay rate usually 
at a level of 5-10T0 -14 . The uncertainty of the phase de- 
lay rate equivalent to group delay, a r /Q„ (where at fl„ is 
the nominal Earth rotation rate) SNR=10 is 300 ps if to 
take into account only thermal noise, and 1400 ps if to con- 
sider the contribution of the atmosphere path delay rate. 
The uncertainty of group delay at SNR=10 in our experi- 
ment is around 200 ps, i.e. a factor 7 better. If to use only 
phase delay rates, the adjustments to nuisance parameters, 
such as atmosphere path delays in zenith direction, station 
positions, and clock functions, will be determined less pre- 
cisely than from the group delay solution and they will affect 
estimate of source coordinates. In order to alleviate the ef- 
fect of poorly determined nuisance parameters on estimates 
of maser positions, I ran a special solution with observa- 
tion equations using both group delay and phase delay rates 
for all detected sources, continuum spectrum objects and 
masers. The phase delay rate had zero weight for continuum 
sources and group delay had zero weight for maser sources. 

Positions of maser sources and their uncertainties are 
presented in Table [2] This table also contains the peak min- 
imum and maximum correlated flux densities for each com- 
ponent of a maser source over observations at different base- 
lines and the velocity with respect to the local system of rest 
(LSR) assuming the barycenter of the Solar system moving 
with the speed of 20 km _1 s towards the direction with right 
ascension 18 h 03 m 50 s .24 and declination +30°00'16".8. The 
table contains the full-width half-maximum (FWHM) of the 
line profile and the minimum and maximum correlated flux 
densities integrated over the line profile: J F(v) dv. For two 
sources that were correlated only with the low spectral res- 
olution, the low limit of their peak correlated flux densities 
and the low limit of the FWHM of their line profiles are 
presented. 

5.1 Comments on individual masers 

Cep A. ICuriel et all (|2002h reported detection of a 2 mjy 
continuum source at 7 mm in the Cep A complex from 
VLA observations in 1996. Position of that source, which 
they labeled as VLA-mm, is 22 h 56 m 17 s .985 +62°01'49".45 
with the uncertainty 0".02. Our position of Cep A water 
maser coincides within l a wit h the position of VLA-mm 
source. iMoscadelli et all l|2009h determined position of the 
methanol maser associated with this object using phase- 
referencing VLBI: 22 h 56 m 18 s . 0970 and 62°0l'49".399 with 
sub-milliarcsecond accuracy. Their position of the methanol 
maser is within 0".78 from the water maser. 

1923+151 is a new discovery. Its spectrum is shown in 
FigureH It is associated with IRAS19230+1506 (3".8 away) 
and WISE J192517.93+151225.0 (0".8 away), within 1-2ct 
of position uncertainty of these catalogues. 

2011+360 is a new discovery. Its spectrum is shown in 
FigureE] It is associated wit h IRAS201 16+3605 (8". 6 away). 
The object was selected bv lSunada et all l|2007f) as a water 
maser candidate, was observed, but not detected. 



6 THE EGAPS CATALOGUE 

Table [3] displays 12 out of 109 rows of the EGaPS cata- 
logue of source positions. The full table is available in the 
electronic attachment. Column 1 contains the class of the 
source, columns 2 and 3 contain the J2000 and B1950 IAU 
names; column 4, and 5 contain right ascensions and de- 
clinations of sources. Columns 6 and 7 contain position er- 
rors in right ascension (without multiplier cos 8) with noise 
4.0 mas/cos S and 1.0 mas added in quadrature to formal 
position uncertainties in right ascension and in declination 
respectively. Column 8 contains correlation coefficients be- 
tween right ascension and declination, column 9 contains the 
total number of observations used in position data analysis. 
Columns 10 and 11 contain the median estimates of the cor- 
related flux density over all experiment of the source in two 
ranges of baseline projection lengths: 0-67 MX (projections 
shorter than 900 km) and 100-170 MX — projections longer 
than 1350 km. If there were no scheduled observations at 
the range of the baseline projection lengths, then —1.000 is 
put in the table cell. 



7 SUMMARY 

I determine positions of 109 compact extragalactic sources 
that are within 6°of the Galactic plane with declinations > 
—20°. Observations were made at a six-element array of the 
Europearn VLBI Network. As a result of these observations, 
the total number of calibrator sources with precisely known 
positions in this region grew from 322 to 431. The semi- 
major axis of position errors range from 5 to 60 mas, and 
the median position uncertainty is 9 mas. 

Among detected sources, nine water masers in our 
Galaxy were identified, two of them new. The date for seven 
these objects were re-correlated with the spectral resolution 
of 0.1 km/s. I determined position of masers from delay rates 
with accuracies 30-200 mas, determined velocities of com- 
ponents and their correlated flux densities. 

I also obtained estimates of source correlated flux den- 
sities. The correlated flux densities of extragalactic sources 
range from 20 to 300 mjy with the median value of 55 mjy. 

The results of this campaign are somewhat disappoint- 
ing. First, the median position accuracy is one order of mag- 
nitude worse than in a similar campaign with the VLBA. 
Several factors contributed to the increase of uncertainties 
in source positions. First, the unfavorable frequency setup 
resulted in an increase of group delay uncertainties by a 
factor of 2.4. Retrospectively, I should admit that more ef- 
forts should be put in order to resolve logistical issues with 
LO frequency changes. Second, the size of the EVN network 
used in this experiment is a factor of 3 less than the size of 
the VLBA network. Third, the sources in EGaPS campaign 
are a factor of 3 weaker than in the VGaPS campaign. 

The detection rate, 30% is significantly lower than in 
other survey experiments. For comparison, the detection 
rate in the VGaPS experiment among sources in the Galac- 
tic plane preselected with VERA observations was 76% and 
36% among the sources selected on the basis of their spectral 
index. Apparently, the efficiency of sources selection based 
on the spectral index is at a level of 30-40% in the Galactic 
plane. About 1/3 detected sources have the correlated flux 
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Table 2. Properties of detected H2O masers. Column description: (1) IAU name, (2) common name, (3) IRAS designator (4) right 
ascension, (5) uncertainty in right ascension in arcsec (without cos<5 factor), (6) declination, (7) uncertainty in declination in arcsec, (8) 
component index, (9) minimum correlated flux density in Jy, (10) maximum correlated flux density in Jy, (11) LSR velocity in km _1 s, 
(12) FWHM of the line profile in km _1 s, (13) minimum integrated flux density in Jy km -1 s, (14) maximum integrated flux density in 
Jy km - 1 s. 
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density below 40 mjy which makes their use as phase cali- 
brator problematic. This indicates that the source selection 
strategy should be revised if a similar survey will be made 
in the future. 
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Table 3. The first 12 rows of the EGaPS catalogue of source positions of 109 target sources. The table columns are explained in the 
text. The full table is available in the electronic attachment. 
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